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ABSTRACT 

An analytical calibration of the Stable Poisson Loaded (SPL) specimen is presented. The specimen 
configuration is similar to the ASTM E-561 compact-tension specimen with displacement controlled 
wedge loading used for R-Curve determination. The crack mouth opening displacements (CMOD’s) are 
produced by the diametral expansion of an axially compressed cylindrical pin located in the wake of a 
machined notch. Due to the unusual loading configuration, a three-dimensional finite element analysis 
was performed with gap elements simulating the contact between the pin and specimen. In this report, 
stress intensity factors, CMOD’s, and crack displacement profiles, are reported for different crack lengths 
and different contacting conditions. It was concluded that the computed stress intensity factor decreases 
sharply with increasing crack length thus making the SPL specimen configuration attractive for fracture 
testing of brittle, high modulus materials. 


INTRODUCTION 

The ability to promote stable crack growth in brittle materials is hindered by the flexibly of the 
loading fixture and high stiffness to fracture toughness ratio of the test specimen. By increasing the crack 
front width with increasing crack length, the Chevron-Notched crack geometries (fig. 1(a)) offer a solu- 
tion to initiating and growing a stable crack in ceramic materials (ref. 1). Unfortunately, the variation in 
the crack width with increasing crack length makes an R-curve type of analysis tedious (ref. 1). An 
alternative technique for establishing stable growth is displacement controlled wedge loading (ref. 2). 
Stability is accomplished by imposing a theoretically constant displacement at the crack mouth causing a 
decrease in stress intensity factors with increasing crack length. The conventional wedge loading configu- 
ration (fig. 1(b)) for a double beam was experimentally attempted and found to be ineffective in produc- 
ing a stable crack in brittle materials with high elastic moduli because the loading train is in series with 
the specimen. The strain energy stored in the loading fixture, which is released during fracture, is often 
enough to cause a catastrophic failure. A solution for minimizing the effect of stored energy is to taper 
the specimen height (fig. 1(c)). The tapered section amplifies the trend of decreasing stress intensity 
factors with increasing crack length (ref. 3). However, the tapered angle chosen must be a function of the 
loading train stiffness (specimen included) and the specimen fracture toughness. In addition to this 



geometric tailoring, side grooves are often used to artificially produce rectilinear crack growth. The exact 
effect of side grooving on measured fracture strength properties is still unknown and debatable (ref. 4). 

A recently introduced displacement controlled wedge loading alternative which produces stable crack 
growth in brittle materials is the Stable Poisson Loaded (SPL) configuration developed by Calomino and 
Brewer (ref. 5). The load train stiffness is isolated from the specimen by loading a cylindrical pin that is 
fitted in the wake of the notched specimen (fig. 2). The specimen is wedge loaded by the diametral 
Poisson expansion of the compressed cylindrical pin. This novel loading configuration has been applied 
successfully to promote stable crack growth in ceramic materials such as aluminum oxide and silicon 
nitride (ref. 6). The specimen geometry is similar to the compact tension geometry and does not require 
side grooving to produce rectilinear crack growth. The loading concept is similar to the ASTM E-561 
compact tension under controlled wedge loading (ref. 7). 

The purpose of this work is to present the numerical calibration of the SPL specimen using the 
finite element method. Two specimen configurations were analyzed: the original SPL with a small 
machined notched width (N » 0), and the modified SPL with a larger notched width (N = 1.77 mm). 

The SPL specimen configuration with N « 0 results in a contacting area equal to 2?rRb, where R is 
the radius of the pin and b is the specimen thickness. Experimental results (ref. 6) indicated that the 
larger pin-to-specimen contact area in the original SPL specimen introduced frictional forces which 
affected the fracture toughness measurements with increasing crack length. To minimize the frictional 
effect at the pin-to-specimen interface, the machined notch width was enlarged in the modified geometry 
(N — 1.77 mm), which reduced the contacting area to xRb/2. Both specimens were analyzed to provide 
an insight in the stability of the SPL geometries. 


NUMERICAL PROCEDURE 

Three-dimensional finite element meshes were generated to analyze the two SPL specimen geome- 
tries. The ceramic specimen and the steel loading pin dimensions used are shown in figure 2. Two types 
of notched geometries were analyzed. In the original configuration, the machined notch width was 
0.5 mm which was taken to be zero for modeling purposes. In the modified SPL geometry the machined 
notch was enlarged (N = 1.77 mm) to minimize frictional effects. The geometrical symmetry permitted 
one quarter of the specimens to be modeled yielding 280 isoparametric 20-noded elements, 1966 nodes for 
the original SPL specimen, and 302 elements with 2196 nodes for the modified geometry. Typical meshes 
for both the modified and original SPL specimens are shown in figures 3(a) and 3(b), respectively. Only 
one element through the specimen thickness was used since the parabolic displacement allowed with this 
type of element usually gives good solutions for both stresses and displacements (ref. 8). The mesh 
density near the crack tip is relatively high in order to accurately capture the stress singularity. At the 
crack tip, quarter-point elements were employed. The central nodes on the edges adjacent to the crack 
front were moved to one-quarter the length of the edge as measured from the crack tip to capture the 
square-root singularity (ref. 9). The ratio of the length of the singularity element over the crack length 
was always kept equal to 0.0097 for all the crack lengths considered. 

The stress analysis was conducted with five different sets of constraints imposed along the interface 
of the two materials due to the uncertainty of the contact condition between the steel pin and the speci- 
men. Nodes along the contacting interface were duplicated to facilitate implementation of the different 
boundary conditions. One set of these nodes was associated with specimen elements, and another set was 
associated with pin elements. 
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The contact between the pin and the specimen was modeled with five different specimen-to-pin inter- 
face boundary conditions: free, fully tied, radial tying, gap element without friction, and gap element 
with friction. Only one coulomb frictional coefficient of /x = 0.40 was considered in this analysis. The 
simplest specimen-to-pin boundary condition (free) linked only the central node along the contacting area 
of the two bodies considered. A fully-tied boundary condition assumed that all contacting nodal degrees 
of freedom for one body were attached to the opposing nodes for the second body. A radial tying bound- 
ary condition, on the other hand, assumed that only the radial degree of freedom for each pair of nodes 
are coupled after transformation to a local cylindrical coordinate system. Another contacting option con- 
sidered in this study was the use of a special gap-link element provided by the general purpose finite 
element program MARC (ref. 8). This element allows separation or contact to occur between the two 
bodies depending on the forces along the link. If the force is tensile, the link provides for separation of 
the two bodies resulting in no transmitted forces. If the force is compressive, the gap closes, and a radial 
force is transmitted. If the friction option is chosen with the gap elements, coulomb type frictional forces 
are generated as the gap closes. The magnitude of the frictional shear force is linearly proportional to the 
normal force, with /j being the coulomb friction coefficient. Thus, material contact produces a frictional 
shear force tangent to the interface between the two bodies. The original SPL specimen required a total 
of 24 gap elements, while the modified geometry used only eight gap elements. The gap elements were 
implemented only on the contacting corner nodes of each 20-node brick element, along the specimen-to- 
pin interface. The mid-nodes purposely were not linked to prevent oscillation of displacements and force 
a linear displacement variation along the element. 

The stress intensity factor (SIF) was determined as a function of crack length from the crack tip 
opening displacement as well as the J-integral approach. The SIF for the free contacting option was also 
calculated using the compliance method for comparison. The plane strain SIF from the crack-tip opening 
displacement is given by the following equation (ref. 9): 
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where t is the length of the singularity element and and V c are the half opening displacements at 
the 1/4 and corner nodes of the singularity element along the crack front, respectively. The variables E c 
and u c denote the specimen Young’s modulus and Poisson’s ratio, respectively. 

The J-integral is calculated from the change in strain energy due to a small perturbation of the 
nodes around the crack tip (ref. 10). A total of 616 nodes are displaced near the crack tip to determine 
the SIF using the following equation: 
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where AU is the change in strain energy, Aa is the change in crack length, b is the specimen thickness, 
and Gj is the mode I energy release rate. 
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RESULTS 


The finite element analysis was conducted for both the steel loading pin and the high elastic modu- 
lus specimen. The elastic properties for both materials are shown in table I. A uniform axial compressive 
pressure was applied to the pin produced a diametral expansion effectively wedging open the crack. The 
SIF and the crack opening displacements (COD) are determined for different normalized crack lengths, 
a/W, ranging from 0.25 to 0.90. The analysis was conducted for both the original SPL specimen with a 
machined notch width nearly equal to zero and the modified specimen with the larger notch (N = 1.77 mm). 


Load Transfer 

The load transferred to the specimen from the pin’s radial expansion was calculated by summing 
nodal reaction forces in the Y-direction (fig. 3(a)) along the symmetry plane of the pin (the Y = 0 plane 
of the pin). The ratio of the load transferred to the specimen, P , to the applied pin load, P a , is shown 
in figure 4 as a function of a/W for each contact boundary concution. The data displayed in figure 4 
was limited to the modified SPL specimen, as this represents the preferred test geometry. All contact 
conditions between the loading pin and the specimen revealed decreasing transferred loads with increasing 
crack length for a given applied compressive pin loading, demonstrating the fracture stability of the speci- 
men. The decreasing trend in transferred load was consistent with increasing specimen compliance. The 
fully-tied contact option always supplied the largest transferred load. 


Crack Opening Profiles 

COD profiles were evaluated for each contact condition in the analysis of the modified SPL speci- 
men. These results are shown in figures 5(a) to 5(e) as a function of the distance measured from the pin’s 
center. The half crack opening displacement V was normalized with respect to E c (the modulus of the 
ceramic specimen), b, and P 8 p, Local deformations in the area of the pin were evaluated for each con- 
tact option considered. 

The large deformations in the pin region, evaluated for the fully-tied case, represent a stretching of 
the material. These deformations were the result of the imposed contact conditions at the interface. 

Such deformations are unrealistic since large tensile forces are exerted to keep the specimen in contact 
with the pin. Deformations near the pin for the radial-tying case (fig. 5(b)) also suggests the presence of 
tensile forces across the pin-to-specimen interface, although the magnitudes are considerably lower than 
that for the fully-tied case. Similar displacement profiles are observed for the same contact conditions in 
the original specimen, although not shown. 

The smoothest crack opening profile displayed (fig. 5(c)) was produced by the free contact case, 
where only the central node of the pin is attached to the specimen. Very similar crack opening profiles 
were seen for the gap element solutions (fig. 5(d)) without friction and (fig. 5(e)) with friction. These 
three contact conditions result in small changes in local deformations near the pin and do not alter the 
overall crack opening displacement profiles. All three crack opening profiles compare well with each 
other. 

Normalized half opening displacements at two specific locations along the crack front are plotted in 
figures 6(a) and 6(b). The displacement, designated Vj in figure 6(a), refers to the displacement at the 
free edge of the specimen located behind the pin at a distance x = — 10 mm. This displacement is also 
referred to as the crack mouth opening displacement (CMOD). The displacement denoted V 2 in figure 6(b) 
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is the crack opening evaluated at a distance x = 7.0 mm (a/W — 0.2) in front of the pin. The fully-tied 
case results in smaller opening displacements than the other contacting options, which is an effect of the 
tensile forces across the contacting bodies. Both displacements, V x and V 2 , increase with increasing 
crack length for all considered contact conditions. The normalized opening displacements of the radial- 
tying and the free case were almost identical for both locations. The opening displacements for the gap 
elements with- and without-friction were slightly lower. They are, however, still higher than the fully 
tied case. Since the load transferred to the specimen is difficult to measure experimentally, the ratio of 
V 1 over V 2 as plotted in figure 7 can be used to obtain the crack length through nonvisual techniques 
as suggested by ASTM E-561 standard (ref. 7). 


Stress Intensity Factor 

The stress intensity factors for the modified specimen with the free contact condition case are given 
in table II as a function of crack length. Listed in table II are the results of three techniques used to 
calculate SIF values: the J-integral method, the crack opening displacement in the singular element 
immediately adjacent to the crack tip, and finally the compliance method. The SIF from the compliance 
method was determined from the change in the work done by the external force with increasing crack 
length. For the free case, the normalized compliance was taken at the only node of the specimen attached 
to the pin. The normalized SIF factor from the compliance method is then given by: 


Y = 


K x b 


sp 



( 3 ) 


where V is half the opening displacement at the center of the contacting bodies, and the derivative of 
the compliance is determined here numerically using the central difference formula (ref. 11). 

As seen in table II, the variation in the normalized SIF between the J-integral and the displacement 
methods is less than 4 percent for a normalized crack length ranging from 0.4 to 0.9. The difference 
increases to 9.7 percent for the shorter crack length. The difference in the SIF, comparing now the 
J-integral approach to the compliance method, is less than 4 percent for a/W ratios between 0.30 
to 0.80. The difference increases to 9.6 percent at the larger crack length, a/W = 0.85. The compliance 
method is affected by the sharp increase in displacement near the free edge. 

In summary, the variation between the numerically determined SIF values from the three different 
methods is within 6 percent for a/W ranging from 0.3 to 0.8, demonstrating a confidence in the results 
within this range. Discrepancies at short and long crack lengths were attributed to the free edge effects 
and the sharp gradient in the opening displacement for the deep cracks. 

Having shown similarity in the results using the three different methods for the free contact option, 
the values of the SIF for all other contact options were calculated for the modified SPL specimen and are 
presented in figure 8. Since the J-integral method has proven to be highly accurate even with a coarse 
mesh (ref. 10), it is the preferred method for computing the SIF results. The reported SIF values were 
normalized with respect to the load transferred to the specimen P gp and evaluated using only the 
J-integral approach. 
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The SIF values for the fully-tied case were always lower than the other contact options. This trend 
was expected since this fully tied case exerted tensile forces between the pin and the specimen. Also 
shown in the figure 8 are the SIF results for the Compact Tension (CT) specimen from reference 12. The 
comparison is adequate since the geometry of the SPL specimen is very similar to the CT specimen. The 
excellent agreement between the present analysis and the CT specimen increases the confidence in the 
numerical procedure adopted here. 

Since the load transferred to the specimen from the pin is experimentally difficult, if not impossible, 
to measure, the SIFs were normalized with respect to the E c and the free edge half opening displacement 
Vj, (Kj/W/E c V 1 ). These are re-plotted in figures 9(a) and 9(b) for the modified and the original speci- 
mens, respectively. As seen from the figures, there is a sharp drop in the stress intensity factor with 
increasing crack length under a constant applied displacement condition. The reduction in stress intens- 
ity factor was about 75 percent, making the SPL specimen geometry very attractive for fracture tough- 
ness testing of brittle materials. The variation in the SIF for the modified specimen was much smaller 
than the original specimen for a given a/W ratio, for the three contacting options shown. When the 
notched width was small (N « 0), the pin introduces additional constraints on the free opening of the 
specimen. This additional constraint is responsible for much wider SIF variations at smaller crack 
lengths. The variation in the SIF for the original specimen between the free and radial tying were iso- 
lated as the source of large experimental scatter in the fracture toughness measurements (ref. 6). The 
modified specimen provided more consistent SIF results, irrespective of the contacting options used as 
seen in figure 9(a). 

The numerical results needed for experimental calibration of the modified SPL specimen are given in 
table III assuming the frictional gap contacting condition with // — 0.4. Although the actual frictional 
forces at the pin-to-specimen interface are not known, its effect is felt to be well modeled since the solu- 
tions are bounded within a narrow region. The gap element option with friction is considered the most 
realistic simulation of actual experimental contacting conditions. Shown in table IV are the best fit equa- 
tions of the results presented in table III for experimental applications. 


CONCLUSIONS 

A three-dimensional stress analysis was conducted on the original and the modified SPL specimens. 
The stress intensity factor calculation for the modified specimen configuration, with a reduced pin-to- 
specimen contact area, showed relatively small dependence on the interaction between the loading pin and 
the specimen. The SIF of the original specimen design, which incorporated a larger contact area between 
the loading pin and the specimen, showed a larger sensitivity to the type of contact assumed. 

By increasing the contact area, the pin introduces additional constraints on the free opening of the 
specimen. This observation is consistent with the large scatter observed in the previously reported frac- 
ture toughness measurements, which utilized the original SPL specimen design. 

The modified SPL specimen configuration is very stable and attractive for fracture testing of brittle 
materials since the SIF is reduced by more than 75 percent for a crack extension from a normalized crack 
length, a/W of 0.25 to 0.90, under constant displacement conditions. 
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TABLE I.— ELASTIC MATERIAL PROPERTIES 
ASSUMED 



Modulus, E 
GPa 

Poisson ratio, 

i/ 

Brittle specimen 

350 

0.25 

Loading pin 

207 

0.30 


TABLE II.— NORMALIZED SIF FOR THE MODIFIED SPL SPECIMEN 
ASSUMING A FREE CONTACTING OPTION 


Normalized 

crack 

length 

Normalized stress intensity factor, 

K,b/w” 

P .P 

a/W 

J-Integral 

Displacement 

Compliance 

0.25 

5.351 

4.833 

-9.7%' 



0.30 

5.905 

5.565 

-5.7% 

5.749 

-2.6%’ 

0.35 

6.579 

6.202 

-5.7% 

6.411 

-2.6% 

0.40 

7.412 

7.181 

-3.7% 

7.233 

-2.4% 

0.45 

8.460 

8.327 

-1.6% 

8.269 

-2.3% 

0.50 

9.807 

9.744 

-0.6% 

9.605 

-2.1% 

0.55 

11.577 

11.565 

-0.1% 

11.372 

-1.8% 

0.60 

13.973 

13.998 

0.2% 

13.777 

-1.4% 

0.65 

17.320 

17.367 

0.3% 

17.350 

-0.2% 

0.70 

22.201 

22.238 

0.2% 

22.190 

-0.1% 

0.75 

29.737 

29.693 

-0.1% 

30.017 

0.9% 

0.80 

42.370 

42.036 

-0.8% 

43.977 

3.8% 

0.85 

66.461 

65.141 

-2.0% 

72.825 

9.6% 

0.90 

123.978 

119.415 

-3.7% 




Indicates percent difference from the J-integral values. 
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TABLE III.— CALIBRATION OF THE MODIFIED SPL SPECIMEN WITH A 
CONTACTING COULOMB FRICTION COEFFICIENT EQUAL TO 0.4 


a 

W 

V, 

V 2 

p .p 

E c V, b 

E c V 3 b 

K, b ffF 

K[ £T 

p. 

P .p 

P -p 

p .p 

E c V 1 

0.25 

4.082 

1.17E-2 

7.98 

1.96 

5.334 

0.6684 

0.30 

3.128 

1.43E-2 

10.12 

3.23 

5.881 

0.5814 

0.35 

2.727 

I.19E-2 

11.49 

4.21 

5.916 

0.5148 

0.40 

2.491 

9.85E-2 

15.93 

6.40 

7.359 

0.4619 

0.45 

2.329 

8.06E-3 

20.03 

8.60 

8.382 

0.4185 

0.50 

2.206 

6.51E-3 

25.38 

11.51 

9.690 

0.3817 

0.55 

2.108 

5.15E-3 

32.62 

15.48 

11.403 

0.3495 

0.60 

2.026 

3.98E-3 

42.81 

21.13 

13.702 

0.3201 

0.65 

1.956 

2.97E-3 

57.82 

29.56 

16.883 

0.2920 

0.70 

1.896 

2.13E-3 

81.23 

42.84 

21.444 

0.2640 

0.75 

1.844 

1.44E-3 

120.35 

65.27 

28.299 

0.2351 

0.80 

1.798 

9.07E-3 

191.84 

106.67 

39.251 

0.2046 

0.85 

1.759 

5.13E-3 

338.92 

192.71 

58.208 

0.1717 

0.90 

1.724 

2.50E-3 

693.26 

402.14 

93.892 

0.1354 
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TABLE IV— CURVE FITS OF THE SIFS AND CODS FOR THE MODIFIED 
SPL SPECIMEN WITH A CONTACTING COULOMB 
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(c) Tapered cantilever beam. 

Figure 1 . — Various fracture specimen configurations for brittle 
materials. 



(b) Original specimen, a/W = 0.5. 


Figure 3.— Finite element mesh for the stable poisson loading 
specimen. 



Normalized half opening displacement, EVb/P 
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(b) V 2 measured at x = 7 mm. 

Figure 6.— Normalized halt opening displacements as a function of 
crack lengths for different contacting options for the modified 
SPL specimen. 



Figure 7. — Ratio of the opening displacement V! over V 2 at two 
different locations as a function of the normalized crack length 
for the modified specimen. 



Figure 8. — Normalized stress intensity factor with respect to the 
load transferred to the specimen as a function of crack lengths 
for different contacting options for the modified SPL specimen. 



(b) Original SPL specimen. 

Figure 9. — Normalized stress intensity factor with respect to the free 
edge opening displacement as a function of crack lengths for 
different contacting options. 


13 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time tor reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503 


1. AGENCY USE ONLY ( Leave blank ) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 


October 1992 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Numerical Calibration of the Stable Poisson Loaded Specimen 


6. AUTHOR(S) 

Louis J. Ghosn, Anthony M. Calomino, and Dave N. Brewer 


WU -505-63-5 B 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-7335 


9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 105609 



11. SUPPLEMENTARY NOTES 

Louis J. Ghosn, Sverdrup Technologies Inc., Lewis Research Center Group, 2001 Aerospace Parkway, Brook Park, Ohio 44142; 
Anthony M. Calomino, Lewis Research Center, Cleveland, Ohio and Dave N. Brewer, Propulsion Directorate, U.S. Army Aviation 
Systems Command, Lewis Research Center, Cleveland, Ohio. Responsible person, Louis Ghosn, (216) 433-3249. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 39 


1 3. ABSTRACT ( Maximum 200 words) 

An analytical calibration of the Stable Poisson Loaded (SPL) specimen is presented. The specimen configuration is 
similar to the ASTM E-561 compact-tension specimen with displacement controlled wedge loading used for R-Curve 
determination. The crack mouth opening displacements (CMOD’s) are produced by the diametral expansion of an 
axially compressed cylindrical pin located in the wake of a machined notch. Due to the unusual loading configura- 
tion, a three-dimensional finite element analysis was performed with gap elements simulating the contact between the 
pin and specimen. In this report, stress intensity factors, CMOD’s, and crack displacement profiles, are reported for 
different crack lengths and different contacting conditions. It was concluded that the computed stress intensity factor 
decreases sharply with increasing crack length thus making the SPL specimen configuration attractive for fracture 
testing of brittle, high modulus materials. 


14. SUBJECT TERMS 


Finite-element; Stress intensity factor; Brittle materials; Fracture toughness testing; 
Stable specimen 


15. NUMBER OF PAGES 


' 14 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Sid. Z39-18 
298 102 























